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Abstract 

We have investigated the influence of different non-linear relativistic mean field models (NL, NLp 
and NLp5) on spallation neutrons for p+Fe and Pb reactions at 0.8, 1.2 and 1.6 GeV by means of a 
relativistic Boltzmann Uehling Uhlenbeck (RBUU) approach plus a statistical multifragmentation 
(SM) decay model. We find that the "evaporation shoulder", i.e. the neutron energy spectrum 
from 3 to 30 MeV, almost for any emission angle is quite sensitive to the isospin part of the mean 
field. For the more neutron-rich Pb-target the evaporation component can be directly related to 
the low density behavior on the symmetry energy in the thermal expansion phase of the excited 
compound system. It turns out that the spallation data for the reactions under study are shown 
to be more consistent with RBUU+SM employing the NLp effective lagrangian. 

PACS numbers: 25.40. Sc, 24.10.Lx, 25.40.Ep, 24.10.Pa 
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I. INTRODUCTION 



Spallation reactions are very important for their wide applications in accelerator technol- 
ogy, production of energy, astrophysics and transmutation of nuclear waste 
The theoretical description of the spallation process is essential in understanding better the 
physical mechanism of each of the mentioned cases. 

Several dynamical models have been constructed for the theoretical description of spal- 
lation process 
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11]. All of them have the same basis, they describe the reaction 



as a cascade of nucleon-nucleon (NN) collisions, but employing different assumptions. The 
main difference concerns implementing mean field dynamics. One can distinguish the sim- 
plest models, which neglect features of the mean field dynamics and employ constant static 
potential, like a class of Intra-nuclear cascade (INC) models [7]. Other, more sophisticated 
approaches comprise changing field and real fluctuations obtained due to use of two and 
three body potentials e.g., quantum molecular dynamics (QMD) models [h]]. Boltzmann- 
Uehling-Uhlenbeck (BUU) type models use both non-relativistic and relativistic mean field 



potentials 0, [l^]. An important advantage of relativistic mean field models is the distinct 
separation of scalar and vector Lorentz-components of the nuclear mean field potential. The 
scalar part is associated with the in-medium dependence of the effective mass (Dirac mass), 
while the vector one changes the properties of the particle momenta inside the hadronic 
environment. This separation between scalar and vector self energies is very useful when 
extending mean field models to asymmetric nuclear matter, in which the role of isovector 
mesons can be clearly associated with the in-medium properties, symmetry potentials and 
effective masses of protons and neutrons. 

In this work, we study the double differential cross sections (DDCS) of emitted neutrons 
in the reactions induced by 0.8, 1.2 and 1.6 GeV proton on targets Fe and Pb by using a 
covariant transport model of a relativistic BUU (RBUU) type 12|. The RBUU model is 
found to be a useful tool in describing the time evolution of the nucleus-nucleus reaction 
dynamics B.QQQ. Asa genuine feature of transport theories it has two in— 
ingredients: a relativistic mean field (RMF) based on quantum hadrodynamics 17J and 
isospin effects that account for the elastic and inelastic channels in resonance production 
and decay (cf. Refs. Q] for different numerical realizations). In our approach, the non- 
linear (NL)-RMF mean-field in RBUU consists of isoscalar and isovector parts with different 
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Lorentz properties. The isoscalar part is characterized by an attractive scalar (a-meson) and 



a repulsive vector fie 
of nuclear matter 



d (cj-meson), which is of importance in describing saturation properties 
17l ] . The isovector part of NL-RMF is also characterized by a competition 
between vector (p-meson) and scalar (#-meson) fields, which is responsible for the density 
dependence of the symmetry energy [ly, [l9|, [20|, |2l| . It should be noted that, the scalar nature 
of the 5- meson leads even to a new interesting effect, important in the reaction dynamics, 
the splitting of the neutron and proton effective masses in isospin asymmetric systems. The 
numerical values of the field parameters, fixed to nuclear matter and symmetry energy at 



saturation, are chosen simi 
results presented in 
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22 



ar to the QMD parameters for a better comparison with the 



23|. 



An important feature of our transport model is that isospin effects are also explicity 



included in the collision integral, see later. T 



RBUU code presented in Refs. P, Q, Q, 



re results presented here are based on the 
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2JJ, where also the properties of the NL- 



RMF, cross sections and the collision integral are discussed. For a meaningful comparison 
with spallation neutron data we use the statistical multifragmenation (SM) model as an 



afterburner 



24]. 



The paper is organized as follows. Section 2 defines the basic ingredients of the RBUU 
model. In Sec. 3, we apply the RBUU+SM code systematically to the experimental neutron 
DDCS in the interactions of 0.8, 1.2 and 1.6 GeV proton on Fe and Pb at fixed angles of 
0°, 10°, 25°, 40°, 55°, 85°, 100°, 115°, 130°, 145° and 160°. We summarize and conclude this 
work in Sec. 4. 



II. DESCRIPTION OF THE RBUU MODEL 



In this section, an outline of the RBUU approach is given, which is described in detail in 



Refs. 
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211 ] . The RBUU equation has the form 



[p n % + (PIF^ + m*«9>*) flf] f(x,p*) = T coll 



(1) 



Eq. [T] describes the evolution of the single particle distribution function f(x,p*) under the 
influence of a mean- field, which enters via effective masses m*, effective momenta p* M and 
the field tensor F^ v = d^Y? — <9"£ M , and 2-body collisions l co u- S s and S M are the scalar 
and vector self energies, respectively. 
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In actual simulations the test particle (relativistic Landau Vlasov (RLV)) method [251.126] 
is used for the numerical treatment of the Vlasov part, the lhs of Eq. ([1]). In this work we 
use a covariant Gaussian (in coordinate and momentum space) shape for the test particles. 



In 25j it was shown that the use of a Gaussian shape for the test particles is appropriate 
to produce smooth fields and it is possible to determine local quantities, such as densities, 
currents, etc, without introducing additional grids. 

With the introduced Gaussians the one-body phase space distribution function f(x,p*) 
is simulated in the following way 

-i A-Ntest °? 

f(x,p')=N—Y, drg (x - x^r)) g (p* - p*(r)) (2) 

— oo 

where A is the number of nucleons and N test is the number of test particles per nucleon. 
The four dimensional Gaussian weights of the test particles take the form 

g(x - Xi(r)) = a s exp (R it ,(x)R^(x)/a 2 ) 5 [(x p - x^(r))wf (r)] (3) 

and 

9{f ~ P*(t)) = * P exp ((p* - P *(T)) 2 /al) 6 [ptf"{ T ) - mf ] (4) 

with a s = (y/TT<7 s )~ 3 and a p = {y/lia p )^. The widths a s and cr p are kept constant, a s is fixed 
normalizing the space Gaussian to unity and a p is correlated making use of the uncertainty 
relation to a s ■ a p = H/2. In Eq. (J3j), 

R?(x) = (x" - *f (r)) - (x„ - x iv {r)) <(r)<(r) (5) 

is the projection of the distance x — Xi(r) on the hyperplane perpendicular to the velocity 
Mj(r). r refers to the eigentime of the particle. 

A non-linear QHD model is adopted for the relativistic mean field potential, with 
isoscalar, scalar a- and vector cu-meson fields and with inclusion of the isovector channel 
through the exchange of the virtual charged , scalar S- and, vector p-mesons 
In this model, in the mean field approximation, the self energies S s and E M are proportional 
to the expectation values of the isoscalar and isovector fields with coupling constants g a , g w , 
g p and g$. The scalar and vector components are given by 

Sf = g^{x)±g p U\x) (6) 

^si = ga(?(x) ±gs5(x) (7) 



19 



2d, 



2i|. 



with the upper (lower) sign corresponds to the proton p (neutron n). The self energies 
characterize the in-medium properties of the nucleons inside the hadronic environment in 
terms of kinetic momenta and effective masses 



v? = v»-n (8) 

m* = M- S si (9) 

which is different for protons and neutrons due to the appearance of the isovector p and 8 
mesons. 

In the local density approximation, the scalar and vector meson fields, determined by the 
scalar density p s and the baryonic current J M , respectively, result from the solution of the 
corresponding equations 

,2„/^ , „„2/ \ , r 3/ n _ „ „ / \ f A A* m *{ X ) 



m a a(x) + acr (x) + bo (x) = g^ps^x) = 9a J d P* r j^\ f( x 'P*) ( 10 ) 

co,(x) = ^J,{x) = ■% / d 4 p%f(x,P*) (11) 
mf, mf , I p 

w = = £j / ^;/3(x, P *) (12) 



(14) 



The isospin vector phase-space and scalar densities are given by / 3 = f p — f n and p s3 = 
p sp - Psn, respectively. 

The equation of motion for the test particles are obtained by substituting Eq. 0) into 
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Eq. (pQ) and putting the collision term equal to zero, such that 
d 



dr 



xt(r) 



U; 



d 



A-Nte 



a s exp(R 2 j (x i )/al) 



- 9a 



± 



m*(xi) 
da(x 



, a t m 
3=1 s 



Z-Ntest 

E 



N, 



test 



m*(xi) * — ' a 
\ l > j= i 



Am 2 



U; 



Am 2 



)) 



a s exp(R 2 j (x i )/al 



Ni 



test 



m*(xi) 



A-Ntes 

E 

j=Z-N+l 



Am 2 



(Rj(xi) - uf UiR jv (xi)) 



a s exp(R 2 (xi) / a\ 



(15) 



test 



Eq. (|T5l) is the Vlasov term of the transport equation in terms of the test particle represention 
of the phase-space distribution function. It describes the propagation of the test particles 
under the influence of the nuclear mean field, which enters here via the coupling functions 

2 

fn — Jjr ( n = Qi 00 an d n — p, 5, for the isoscalar and isovector parts of the nuclear potential, 
respectively). For charged baryons the Coulomb force represents an additional term in the 
equation of motion. 

In the RBUU code, not only protons and neutrons are propagating separately according 
to their hadronic fields and Coulomb interaction but also A 0,± ' ++ -resonances. iV*-resonances 
are not fully accounted for, which may set the upper limit of the incident energy of the code 
to 1.5 GeV for the nucleon induced reactions. The self energies of the resonances are built 



as 
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22 1 



Si(A-) 
S,(A°) 

S,(A+) 
S,(A ++ ) 



Ej(n) 

2 1 

-S«(n) + -£i(p) 

1 2 

-E;(n) + -E;(p) 

£i(p) , 



(16) 



where i = scalar, vector . 
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Differently from A and nucleons, pions are propagated as free particles with respect to 
the nuclear mean field, but experience the Coulomb potential generated by all the other 
charged particles. 

The collision integra 



Carlo algorithm 
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i.e. the r.h.s of Eq. ([T]) is modeled by a parallel ensemble Monte 
16l . |28I |. A geometrical interpretation of the cross section is used, 
so that in each ensemble two test-particles are allowed to collide if their relative distance d 
fulfills the relation: 

d< \I^F L ~ ■ ( 17 ) 

The collision is accepted if the final state Pauli blocking is not violated. 

The total cross section a tot depends on the center of mass energy and on the species of 
the incoming particles. The distance d is calculated using the covariant calculation proposed 



by Kodama et al. 30]. 



The RBUU code describes the propagation and mutual interactions of nucleons, Delta 
resonances and also 7r-mesons. The following hadronic reactions are treated explicitly: 

• NN —>■ NN (elastic scattering) 

• AA — > AA (elastic scattering) 

• NN <-> iVA 

• Nti <-» A 



The different NN «-> NA and iViV <-> A A isospin channels are treated according to [3J|, 
in which the cross sections are evaluated on the basis of the One Boson Exchange model 
and then fits to the analytical calculations are given. For the elastic cross sections and the 



331 ] is used. 



angular distributions, the free parametrization according to Cugnon et al. 
Besides production and absorption in iVW-channels, resonance states can be populated 
also through 7rN collisions. It should be mentioned that, in contrast to other transport 
approaches , [loj], the RBUU code accounts for both the relativistic kinematics and 
the isospin effects in resonance production and decay, due to the Lorentz structure of the 
isovector channel. 

The mean field is determined by four coupling constants, /, = gf/mj, i = a, u, p, 5 
and the two parameters of the non-linear self-interactions of the cr-meson field a, b. In this 
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L (fm~ 2 ) 


L (fm- 2 ) 


f P (fm~ 2 ) 


fs (fm- 2 ) 


B (fm' 1 ) 


C 


NL 


9.3 


3.6 


0.0 


0.0 


0.015 


-0.004 


NLp 


9.3 


3.6 


1.22 


0.0 


0.015 


-0.004 


NLp5 


9.3 


3.6 


3.4 


2.4 


0.015 


-0.004 



TABLE I: Coupling constants fi 



for i = a, u, p, 5, B 



4 and C 



\ for the different 



NL - RMF models. 

work, for the isoscalar part the parameter set NL of Ref. [32] has been used. For symmetric 
matter it gives reasonable values for the incompressibility K = 200 MeV with a nucleon 
effective mass of m* = 0.83M at saturation density po — 0.148 fm~ 3 . In order to unravel 
the role played by the different density behavior of the symmetry energy, we will discuss the 
results obtained with the two parameter sets NLp and NLpS in the isovector channel. In 
NLp and NLpS models the symmetry energy parameter is fitted at saturation to the value 
of 30.7 MeV. The values of all the coupling constants are given in Table 1. In Fig. [1] we 
plot the density dependence of the symmetry energy for the NL, NLp and NLp5 models. 
We note two features that will be relevant for spallation reactions: i) The symmetry energy 
is in general larger when isovector bosons are inserted (in the NL model only the kinetic 
Fermi contribution is present); ii) The contribution of the 5-meson has different effects on 



the symmetry term, a decrease below saturation and an increase above po 3, llSl, 120 ] . 



It should be noted that the NL, NLp and NLpS parameters have been already used 



for pion/kaon production and isospin tracer ca 



with an overall agreement with the data 
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culations in relativistic heavy ion collisions 



28 



29|. 



Event Characterization 



Due to the stochastic nature of the RBUU simulations the transport code is used as an 
event generator. In order to reduce statistical fluctuations 500 events are analyzed for each 
initial condition (target, beam energy and centrality). 

An important issue for the description of spallation data within the RBUU model is an 
appropriate initialization of the ground state nuclei. We will follow here in particular the 
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initialization of Ref. 



31] , which aims to construct the static properties of the ground nu- 



clear state such as binding energy and root mean square radius. Moreover the evolution of 
these quantities should be stable enough, and the behavior of the density and momentum 
distributions should be approximately correct. To get the initial nuclear phase space dis- 
tribution, we first distribute the N test x A test-particles according to a Saxon- Wood shape 
with neutron (proton) radius R n ^(A) and diffuseness parameter a ntP (A) depending on the 
target mass A. The parameterizations of R n ^ p (A) and a niP (A) as a function of mass A for 
stable nuclei have been deduced from the adjustment of the Saxon- Wood function on the 
density matter distributions calculated with the model of Ref. [13J. The Saxon- Wood tail 
is cut off at Rmax = R n , P + 2a„ iP . The next step is to determine the local density of all 
nucleons generated by all the other nucleons. The initial momenta of N test x A test-particles 
are randomly chosen between zero and the local Fermi momentum: p™ ax = hc(3n 2 psax) , 
with psax being the corresponding Saxon- Wood -neutron or proton density. Finally, an ini- 
tial ground state configuration is accepted when (a) the sum of the total energy is equal 
to Eb/A ± 0.5MeV/A, where Eb/A is the ground state energy of a nucleus of mass A and 
charge Z given by the liquid drop model and (b) the binding energy and root mean square 
radius keep smooth variation with time and without spurious particle emission. 

The RBUU calculation is carried out up to a time scale referred to as the transition 
time ttr, when an excited compound system is formed, the "prefragment" source that will 
decay via statistical emissions. At tt r , the position of each nucleon is used to calculate the 
distribution of mass and charge. The minimum spanning tree method [35] is employed, a 
prefragment is formed if the centroid distances are lower than R c ius- In this paper, R c ius 
is fixed at 4 fm. The prefragment thus identified is then Lorentz boosted into his rest 
frame to evaluate the excitation energies. The excitation energy (e exc ) is calculated as the 
difference between the binding energy of the hot prefragment and the binding energy of this 
prefragment in the ground state, evaluated as indicated before. When the prefragment is in 
the excited state, the SM model is employed as an afterburner, which is shown to be 



well suited for the description of the slow evaporated neutrons III, |23|, \34\. 

In order to determine t tr , it must be verified, whether information obtained from the 
transport code are sensitive to the time duration of the first stage of the reaction. For 
this purpose, the time variation of the average values (over 100 events) of three physical 
quantities: excitation energy per nucleon < e exc >, angular momentum < L > and mass 
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number of the excited prefragments < A exc > after the first RBUU stage, have been analyzed. 
As an example, Fig. [2] shows the time evolution of these quantities for p + Pb collisions at 
1.2 GeV, obtained with a NLp interaction. As one can see, large prefragments are produced 
at earlier times 5 < t tr < 10 fm/c with a maximal value of < e exc > MeV/A and relatively 
large < L >. Then, both quantities < e exc > and < L > drop very rapidly in the time 
range 10-35 fm/c, in correspondence of a fast, pre-equilibrium, nucleon emission as shown 
by the decrease of < A exc > in the bottom panel. The dropping of < e exc >, < L > and 
< A exc > stops at around 40 fm/c and the average values remain stable. At later times, 
t > 80 fm/c, spurious increase (decrease) of < e exc > (< L > and < A exc >) is observed, 
which is unphysical and arises due to spurious numerical fluctuations. This may indicate 
that the first dynamical stage of the reaction evolution should be terminated at about 80 
fm/c and a statistical decay procedure could be inserted to get the final reaction products. 

Let us now check the dependence on the t tr choice of the final results. In Fig. [3] we 
show the neutron double differential cross section for the reaction p + Pb at 1.2 GeV. The 
line histograms denote the RBUU+SM calculations using NLp parameterizations with three 
different transition times, 60 fm/c (thin lines), 75 fm/c (thick lines), and 100 fm/c (dot- 
dashed lines). As one can see, calculations with different transition times yield the same 
results in the low energy neutron region rs E n < 30 MeV. This may indicate that the decay 
processes of the excited prefragments described by RBUU+SM are equivalent in a t tr choice 
from 60 to 100 fm/c. This is important for the study of the symmetry term dependence of 
this component of the neutron spectrum presented in the next section. 

On the other hand, in the "cascade region" (see later), 30 < E n < 60 MeV, a systematic 
difference between the calculations is seen as the angle increases. This difference decreases 
within the time interval from 75 to 100 fm/c. We observe that the neutron data, in this 
energy range, are better reproduced by the RBUU+SM at t tr = 75 fm/c. We, thus, use 
this value for all systems in the present study, which is also supported by Fig. [2l 

It should be pointed out here that in the SM model, standard symmetry coefficients 
Isym = 23 to 25 MeV of the fragment binding energy are used to describe light and heavy 
fragments produced in proton/nucleus- nucleus collisions [241 ] . On the other hand, the sym- 
metry energies within the different parameterizations used in the transport model may differ 
from those used in the statistical approach for densities beyond saturation. Therefore, it is 
worthwhile to examine the consequences of changing 7 S3/m in SM calculations, as has been 
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done in different studies that focused on isotopic distributions [3j|-|39j]. Fig. @] displays the 
evaporation neutron spectra obtained using RBUU (with NL)+SM calculations for p+Pb 
reactions at 0.8 GeV, and the symmetry term 7 SJ/m was varied between 5 to 25 MeV. As one 
can see, with 7 SJ/m = 5 MeV, the SM processes causes a slight increase of low energy neutrons 
at E n < 10 MeV. For larger values of r j S ym, however, the evaporation neutron spectra are 
not much affected by the symmetry coefficient difference in the SM code. We note that at 
variance with the Heavy Ion Collisions studies of Refs. [36-38] in our case we do not have 
large compression-expansion effects. The final equilibrated source is then at densities close 
to the saturation value. This has been checked analyzing the time evolution of the root mean 
square radius of the composite system in the (p+Pb) case at 0.8 GeV. After the thermal 
expansion the system shows a rather good stability with smooth monopole oscillations, see 



also Ref. 40j. This may imply that the equilibrating source is formed at about normal 
nuclear density and it is not necessary to use much reduced symmetry coefficients. Thus, 
in what follows the value of 7 s2/m = 25 MeV will be kept fixed in the SM code, which is 
supported by the experimental analysis of ^sym with the SM model in p, d and a induced 



reactions at relativistic energies up to 15 GeV [39 ]. 



III. RESULTS AND DISCUSSIONS 

In this section we present and discuss the predictions of the RBUU+SM model using 
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421 ] of energy-angle DDCS of 



different NL parameter sets along with the measurements 
neutrons induced by 0.8 GeV, 1.2 GeV and 1.6 GeV protons on Fe and Pb. 

The experimental energy spectra (see Figs. [5] and [6] and similar ones) show at 0° two 
prominent peaks. These peaks are less pronounced at 10° and are insignificant at 25° and 
larger. The (quasi-elastic) peak, characterized by a narrow peak at a kinetic energy near that 
of the beam energy, is due to a single (p, n) elastic scattering in the forward direction. The 
(quasi-inelastic) peak, centered around 873 MeV and 760 MeV at 0° and 10°, respectively, 
is about 400 MeV wide and is thought to be due to A-resonance excitation. In addition to 
these two peaks, it seems that two components exist for all of the spectra: one is a "shoulder" 
below E n pa 20 (30) MeV for Fe (Pb), the other is a wide peak extending up to a few hundred 
MeV. The low energy neutrons mostly come from evaporation of the excited target residues 
formed through the equilibration process. The other component becomes less pronounced 
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with increasing angles. This component arises from "cascade" processes involving several 
iVTV-collisions. Below we are going to investigate the effect of the NL, NLp and NLpS 
models on the three components of the neutron spectra by employing RBUU+SM code. 
We have performed 500 simulations at various impact parameters, from to 4(7.5) fm for 
Fe(Pb), respectively. In order to have sufficient statistics, calculations were done for angular 
bins of ±3.5 at 0° and ±5 for larger angles. 

Before confronting RBUU+SM with experimental data, it is worthwhile to investigate 
the sensitivity of the calculations to the number of test particles (N test ) . This has been done 
(not shown here) for p+Fe reactions. We used NLp parameter set and the calculations were 
performed with 40 and 100 test particles per nucleon. It was found that the overall behavior 
of the neutron DDCS are nearly the same for both calculations at all angles. Since the 
computational time of the collision integral scales as (N test x A) (parallel ensemble method), 
a value of N test = 40 test particles per nucleon has been chosen. 

All the results are presented and compared to data in the Figs. loTfTOl Before passing 
to a more detailed analysis we like to note that all the components of the neutron spectra 
are rather well reproduced by this relativistic transport approach. We remark that no free 
parameters are used while genuine relativistic effects are present via the full covariance of 
the model and the natural presence of Lorentz forces due to the coupling to vector fields. 

Let us first focus on the role of symmetry energy on the neutron spectra. This is illus- 
trated in Figs. [MEl by comparing the RBUU calculations without (Figs. [5] and [6]) and with 
(Figs. [7] and E]) evaporation using NL (dot-dashed lines) and NLp (thin lines) models. The 
two models differ only in the isospin part, with p-coupling giving an almost linear density 
dependence of the symmetry energy with value ~ 30.7 MeV at saturation (see Fig. [I]). 

Figs. [5] and [6] show that the results of both calculations, NL and NLp, are close together 
and can describe the high energy part of the neutron spectra (above E n pa 20 and 30 MeV 
for p+Fe and Pb, respectively). It is interesting to see that without statistical decay the 
"low energy shoulder" is absolutely missed. For the more neutron- rich Pb target some 
difference originating from the p-meson becomes more pronounced in the lower part of the 
neutron energy spectrum (below 30 MeV), especially with increasing emission angle. Such 
behaviour can be related to the difference in the symmetry energy: neutrons experience a 
stronger repulsive field going from NL to NLp, at any density (see Fig. [1]) [if], 0]. 

In fact the largest isospin effects are observed just on the low energy evaporation compo- 
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nent which appears when the afterburner statistical decay is included. In Figs. [7] and [8] we 
compare the full RBUU+SM calculations (with evaporation) along with the experimental 
neutron spectra. As one can see, the NL results, with no (potential) symmetry energy, show 
a broad maximum at E n w 3 — 50 MeV in all angular intervals for all studied collisions. 
On the other hand, the NLp calculations reduce the neutron yield in this energy region and 
more closely reproduce the experimental data, especially as the incident energy increases. 
The importance of the isovector interaction is shown to increase when both the incident 
energy and the isospin asymmetry of the target system increase. It should be pointed out 
here that, the effect of the symmetry potential on the three components of the neutron 
spectra has already been discussed by a non-relativistic QMD approach for p+Pb at 1.2 and 



1.6 GeV 



23j. It was shown that of all the three components of the neutron spectra, the 



evaporation part is the most sensitive to the symmetry potential. 

Next we compare, in Figs. [9] and [10J the RBUU+SM results using NLp without (solid 
lines) and with (dot-dashed lines) the 5-field. In the case of p+Fe interactions, there is almost 
no significance differences between NLp and NLp5 results. However, for p+Pb interactions, 
with increasing neutron to proton asymmetry, NLpS generates an excess of slow neutrons 
(E n w 3 - 50MeU) compared to NLp. 

The origin of the different behaviour between the results in the slow neutron energy 
region (j50 MeV) is clearly due to the excitation energy and isospin content of the hot 
prefragments. The properties of these prefragments at the final stage of the interactions 
(at 75 fm/c) are affected by the number of pre-equilibrium neutrons (N em ) emitted at the 
earlier times: less N em leads to larger excitation energy and larger N/Z of the prefragment. 
The NLp calculations of the p+Pb interactions at 1.2 GeV show, in Fig. [5J that nearly 20 
nucleons are emitted at the earlier times (between 20 and 40 fm/c). In Fig.[TT]we examine the 
distributions of emitted neutrons N em before 75 fm/c by employing RBUU with NL, NLp 
and NLpS cases for the p+Pb reactions at 0.8 GeV. One observes that for N em > 8, less pre- 
equilibrium neutrons are emitted with NL < NLpS < NLp. Consequently, the excitation 
energy and isospin asymmetry of the prefragments should appear larger for NL calculations 
with respect to NLp (and > NLpS) as observed. The sequence (NL < NLpS < NLp) 
agrees with the sequence of the three forms of the symmetry term at the densities below 
saturation (see Fig. [[]), which confirms that N em is generated at earlier times where an 
expansion mode, due to thermal pressure, is likely occurred. 
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The sensitivity of the evaporation component of the neutron spectrum to the low density 
behavior of the symmetry term can be then used to disentangle among different effective 
fields and to get an independent information on the symmetry energy below saturation, see 



some general discussions in the reviews 



IS 



431 ] . The best agreement with the data in the 



more n-rich Pb case is obtained within the NLp choice, which corresponds to a rather stiff 
increase of the symmetry energy below saturation, roughly proportional to the total density. 
It is also interesting to notice (see Figs. IBl and [TO]) that the NL and NLpS give similar extra 
yields with respect to the NLp case, which implies that the excitation energies and N/Z of 
the Pb hot prefragments are similar. This can be also explained by considering the weaker 
symmetry term at densities below saturation (see Fig. [T]) 

Finally, let us identify the effects of the mean field on the high energy part of the neutron 
spectra (E n > 40 MeV) for the reactions under study. This is illustrated in Figs. [5] and [6] and 
FigUl] which show a close up of this region at very forward angles (9 < 10°). Two different 
scenarios are explored: RBUU calculations with NLp (thin lines) are contrasted by RBUU 
with deactivated mean field (CASCADE) (thick lines). We observe that at 10° < 9 < 55° the 
NLp calculations are practically identical to CASCADE calculations as the incident energy 
increases, showing the mean field effects to become negligible. At 55° < 9 < 160° the NLp 
calculations give lower values compared to CASCADE calculations, in good agreement with 
the data. However, the NLp calculations underestimate the high energy neutrons (70-110 
MeV) at 9 = 160° for p+Fe at 1.6 GeV and p+Pb at 1.2 GeV, contrary to CASCADE 
calculations. Note that, both calculations resemble each other at 9 = 160° for p+Pb at 1.6 
GeV. Thus, it seems that the high energy neutron spectra evaluated by RBUU (with and 
without NLp) have different beam energy dependence compared to the experimental data. 
It should be pointed out that our results in the backward direction are nevertheless better 

n n 

with our present RBUU code than with the non-relativistic INC [8j, QMD [22| and BUU 
44j calculations. 

Fig. [T2l shows how the quasi elastic and quasi inelastic peaks are changed due to the mean 
field. In the CASCADE calculations, the maximum of the calculated quasi-elastic peak is 
shifted toward a higher energy in comparison to the data. The shape and yield of the peaks 
are in better agreement with the CASCADE calculations at 1600 MeV (and at 10°) than at 
800 and 1200 MeV. On the other hand, as the mean field calculations are taken into account 
in RBUU, the maximum outgoing neutron energy is strongly attenuated at incident energy 
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of 1200 MeV and 1600 MeV. As for the quasi- inelastic peak, i.e., the peak located at the 
beam energy minus ~ 300 MeV, the evolution of the shape with the incident energy follows 
more or less the data for both calculations. We notice that the NLp calculations have a 
tendency to become smaller compared to the CASCADE ones at lower incident energy and 
need to be improved. 

In [45j the quasi-elastic and inelastic peaks at 1.2 GeV were underestimated by the 
ultrarelativistic QMD (UrQMD) calculations, which was attributed to the free angular dis- 
tributions of the scattered NN-elastic and inelastic collisions adopted in the UrQMD model, 
and it was suggested to introduce medium modified angular distributions (MMAD) into the 
UrQMD model for reducing the deviation. We then further test the influence of the MMAD 



on the reactions under study. We introduce the MMAD adopted in [45| in the collision 
term of RBUU (with NLp) code. The dot-dashed lines in Fig{12] illustrate these calcula- 
tions for p+Pb at 800 MeV. As one can see, the introduction of MMAD in RBUU improves 
the intensity of both the quasi-elastic and inelastic peaks. For the location, however, the 
maximum of the calculated quasi-elastic peak is 50 MeV downward shift in comparison to 
the data. At higher incident energy (not shown here), this shift is increased even further to 
100 and 150 MeV for the reactions under study at 1200 MeV and 1600 MeV, respectively. 
We attribute this failure to the linear energy dependence of the Schrodinger equivalent op- 



461 ] . Indeed, Dirac phenomenology 



tical potential of the NL-RMF models at high energies 
on elastic proton-nucleus scattering predicts a nonlinear energy dependence of Schrodinger 
equivalent optical potential for incident energy above 400 MeV, which leads to a softening of 
the optical potential and, consequently, may improve the dynamics of NN-elastic collisions 
at high energies. 



IV. SUMMARY AND CONCLUSIONS 



We have studied proton-induced Fe and Pb reactions at 0.8, 1.2 and 1.6 GeV by a hybrid 
model, using the RBUU transport approach for the first pre-equilibrium stage and the SM 
decay model for the second slow stage of the reaction. The NL-RMF model of the present 
approach consists of isoscalar (a — uj mesons) and isovector (p — 5 mesons) parts with 
different Lorentz covariant properties. The impact of different NL-RMF models (NL, NLp 
and NLpS) on the neutron spectra of the studied reactions has been investigated and the 



16 



following conclusions can be drawn: 

• Some effects of the mean fields are found on the high energy part of the neutron 
spectra, when switching off the mean field and turning to the CASCADE mode within 
the RBUU calculations. 

• The low energy neutron spectrum ( E n pa 3 — 50 MeV) is found to be quite sensitive 
to the details of the relativistic structure of the isovector part of the mean field: the 
inclusion of the p-field improves the comparison between theory and experiment to 
large extent. The effect is related to the properties (rate and isospin content) of the 
fast nucleon emission during the thermal expansion of the system and so it is sensitive 
to the low density behaviour of the symmetry term. 

• The competition between a scalar attractive (5) and a vector repulsive (p) fields affects 
low energy neutrons (E n < 30 MeV) only for the more neutron-rich Pb target. 

• The introduction of medium modified angular distribution in RBUU (with NLp) im- 
proves the intensity of p+Fe and Pb reactions at 0.8 GeV. 

Thus, the present comparisons suggest a better agreement with NLp model, however, 
several other aspects of the mean field dynamics of the RBUU model need further work. For 
example, the nucleon self-energies of the current version of the model have no momentum 
dependencies, which might improve the extremity of the calculated spectra at 0° and 160°. 
This work is in progress. 
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FIG. 1: (color online) Density dependence of the symmetry energy for the NL-RMF models. 
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FIG. 2: (color online) Calculated average excitation energy (< e exc >), angular momentum (< L >) 
and mass number of prefragments < A exc > as a function of transition time (tt r ) after the RBUU 
(with NLp) initiated by a proton on Pb at 1.2 GeV. 
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FIG. 3: (color online) Neutron energy-angle double differential cross sections for p+Pb at 1.2 GeV 
coated by RBUUfcj* N L p)+ SM wita differ tuition ta. M . Data (soM c W e S with 
error bars) are from [41|]. For clarity, only the histograms and the data for the smallest angles are 
given in absolute value. The other ones have been multiplied by 10 -1 , 10 -2 , . . . for other angles 
in increasing order. 
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FIG. 4: (color online) Evaporation neutron double differential cross sections for p+Pb reactions at 
0.8 GeV calculated by RBUU(with NL)+SM. The lines are the results obtained by changing the 
symmetry coefficients {^ S ym) in the SM code. 
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FIG. 5: (color online) Neutron energy-angle double differential cross sections for p+Fe interactions 
at 0.8 GeV(left panel), 1.2 GeV(middle panel) and 1.6 GeV (right panel) as compared to RBUU 
calculations (without evaporation). The dashed and thin lines denote calculations with NL and 



4JJ (solid 



NLp, respectively. The thick lines show calculations without mean field. The data 
circles with error bars) and the histograms for the smallest angles are given in absolute value. The 
other ones have been multiplied by 10 _1 , 10 -2 , . . . for other angles in increasing order. 
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FIG. 6: (color online) Same as Fig. 5, but for p+Pb reactions. Solid triangles with error bars (left 
panel) are data taken from [42]. 
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FIG. 7: (color online) Same as Fig. 5, but the solid and dashed lines denote the RBUU+SM with 
NLp and NL, respectively. 
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FIG. 8: (color online) Same as Fig. 7, but for p+Pb reactions. Solid triangles with error bars (left 
panel) are data taken from [42]. 
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FIG. 9: (color online) Same as Fig. 5, but the dashed and solid lines denote the RBUU+SM results 
with NLp5 and NLp, respectively. 
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FIG. 10: (color online) Same as Fig. 9, but for p+Pb reactions. Solid triangles with error bars (left 
panel) are data taken from [42]. 
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FIG. 11: (color online) Normalized distribution of pre-equilibrium neutrons for p+Pb reactions at 
0.8 GeV. The lines represent the RBUU (without SM) with different NL-RMF models. 
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FIG. 12: (color online) Neutron energy spectrum at 0° and 10° from 0.8 GeV (left panel), 1.2 GeV 
(middle panel) and 1.6 GeV (right panel) proton interactions with Fe (top panels) and Pb (bottom 
panels) targets. The thin lines denote RBUU calculations with NLp, while the thick lines are 
those without mean field. The dashed lines denote RBUU (with NLp) calculations that include 



the medium modified angular distributions of Ref. 



45j. 
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